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ABSTRACT 


A digital computer program for the geometrically 
nonlinear analysis of totally arbitrarily loaded 
shells of revolution (SATANS-II) was modified to more 
accurately account for the conditions at the pole of 
the shell. This program was used to determine the 
buckling load of shallow spherical sheils of various 
sizes when subjected to static axisymmetric, dynamic 
axisymmetric, and dynamic nearly axisymmetric 
step-pressure loads of infinite duration. A 
comparison was made between the new buckling results 
and previous results obtained without the new pole 
routine. The comparison, revealed a significant change 
in the buckling pressures, due solely to the change in 
the pole routine. The new static axisymmetric, 
dynamic axisymmetric, and even the dynamic asymmetric 
critical buckling pressure loads appear to be fairly 
reliable results for perfect, shallow shells. 
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I. INTRODUCTION 


In 1973 a digital computer study was presented by Ball 
and Burt [{1] for the dynamic buckling load of clamped 
shallow spherical shells subjected to axisymmetric and 
nearly axisymmetric step-pressure loads. A static buckling 
analysis of the same spherical shells had been. carried out 
ap 1970 by Stilwell and Ball [2]. In these two studies the 
digital computer program SATANS-I [3] was used to calculate 
the critical buckling pressures for a large range of shell 
sizes. Other studies of the buckling of shallow shells have 
been conducted by Huang {4,5]J, by Stephens and Fulton [6], 
by Lock et al. [7], by Stricklin [8], and most recently by 
Akkas [(9]. In Reference 1 the results from these other 
studies, except for those by Akkas, are compared with the 
results from SATANS-I for both static and dynamic buckling. 
In the axisymmetric static analysis tae comparison with the 
results obtained by Huang [4] revealed that the SATANS-I 
results were higher than Huang's results for several shell 
sizes. In the dynamic, axisymmetric buckling analysis the 
SATANS-I results again either agreed closely with, or were 
somewhat higher than , the results by Huang [5], Stephens 
ana Fulton [6], and Stricklin [8]. However, it was noted 
then that there was a general lack of consistent agreement 
among any of the sets of results. AS a consequence, it 
appeared at that time that the axisymmetric buckling problem 
had not yet been totally resolved and that additional 


studies would be appropriate. 


In the asymmetric dynamic buckling analysis of Reference 
1 the few comparisons that could be made for the critical 
load also indicated that the SATANS-I results may be too 
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hagh. A ccmparison of the- recent estimates for the 
asymmetric dynamic buckling load obtained by Akkas [9] with 
the SATANS-I results also reveals the SATANS-I resuits to be 
well above those of Akkas [9]. However, it should be noted 
that the results obtained by Akkas were from his attempt to 
obtain a lower bound on the critical asymmetric load. This 
bound on the buckling load is obtained without the execution 
of a complete transient response analysis on the asymmetric 
part of the response of the shell, as is done in SATANS-I. 
In Akkas' analysis (Problem 1) the transient nonlinear 
axisymmetric response is computed, and a determinant is 
examined for possible bifurcation into asymmetric moticn at 
each time step. The minimum load at which the determinant 
becomes zero is defined as the lower bound of the critical 
load. 


AS a consequence of the generally high buckling loads 
predicted by SATANS-I, a re-examination of the static and 
dynamic buckling of the shallow spherical shell was made in 
an attempt to determine the possible cause, or causes, of 
the high buckiing loads. In our search we discovered that a 
modification of the manner in which the pole conditions are 
numerically approximated significantly lowered the buckling 
loads to values that are now in good agreement with the 
other results. The new procedure for handiing the pole 
condition is given in section III of this thesis. The new 


buckling results are given in section V. 


In addition to the pole condition modifications and the 
new buckling results the author has also made another 
Significant change to the SATANS family of codes. [In 
particular, the SATANS-II program for the geometrically 
nonlinear analysis of totally arbitrarily loaded shells of 
revolution, developed by Ryan [10] in 1972 to handle more 
complex and larger problems, was modified to make the 
computer memory requirement a variable quantity. This 
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@uautity iS specified by the user to fit the particular 
problem being run. It eliminates the large core requirement 
of SATANS-II for small problems and allows for much larger 
problems to be soived than could be solved by SATANS-II. 
The new program with the pole condition and memory 
modifications will hereafter be called SATANS-IIA. Tee is 


described in section II. 
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J od ESCRIPTION OF SATANS-IIA 


SATANS-II was developed by Ryan [10] from SATANS-I and 
incorporated the full trigonometric expansion of the applied 
load and solution vector, and introduced the handling of 
imperfections into the code. These modifications allow the 
analysis of shells under totally arbitrary loads, as well as 
imperfection studies on actual shells with measured 
imperfections [11]. Unfortunately, the original deck of 
cards for SATANS-II was destroyed. Professor Johann Arbocz 
of CALTECH had a listing of SATANS-II and punched a deck of 
cards with the changes to SATANS-I given in that listing. A 
copy of this deck was sent to Professor Ball. These cards 
have been added by the author to the original SATANS-I 
described by Ryan [10] and a complete version of SATANS-II 
has been reconstructed. SATANS-IIA is a modification by the 
author of the reconstructed SATANS-II program. A listing of 
SATANS-IIA can be found in Appendix A. The listing contains 
an example problem for the dynamic analysis of a clamped, 
truncated cone subjected to an impulsive loading which is 
uniform along the meridian and varies in a cosine 
distribution over one-half of the circumference. This 
problem is a sample problem suggested by the Lockeed 
Missiles and Space Corp. [12]. A condensed version of the 
output from the example problem is given in Appendix B. 
Input data preparation for SATANS-IIA can be found in 
Appendix C. The basic users manual, which includes 
preparation of input subroutines and the theory of the 
program, is contained in Reference 3, which can se obtained 
toEough COSHNIC..(470-10098,. LAR-10736), or ASIAC [13]. A 
users manual which includes preparation and handling of 
imperfection data within the SATANS programs can be found in 
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Ref. {10]. The above information, along with the following 
discussion, will inform the user on the capabilities and 


proper use of SATANS-IIA. 


The modification of the SATANS-II program to make its 
core requirement variable was accomplished by putting in a 
single dimension statement at the beginning of the progran, 
with subsequent dimensioning within the subroutines to only 
ne first “element -of-~the vector or matrix. This is a 
convenient feature of the FORTRAN-IV language in which the 
program is written. The actual vector and matrix sizes are 
transmitted to the subroutines by. an individual parameter 
EESt. Construction of the initial dimension statement and 


core request size is as follows: 


The basic size of the program on the IBM-360/67 Digital 
Computer, without the initial dimension statement, is 
272,000 bytes. This figure includes approximately 19,000 
bytes of buffer space required for execution. Within the 
Main dimension statement are fifteen variables. However, 
only three parameters are needed to specify the sizes of 


these fifteen variables. 


Let a= The number of stations along the meridian of the 
shell times the number of harmonics considered. 


Let b= a, plus two fictitious stations times the number 


of harmonics considered. 
Let c= The number of harmonics considered. 


The main dimension statement would then be constructed 
as, 
DEMENSTON P(l Goa), DER(4, 4, a), DST(4,4,a), X(4-a) » 
PHEXB(a), PHITB (a), Z2(4,b), 2ZO0(4,b), 
Z2teb),) Star b) , ZDOT (4,b)» 1S9(99,c) , 
JS(99,c), 1D(99,c), JD(99,Cc) 
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The 99's above limit the user to 99 karmonics in any one 
run and an unlimited number of meridional stations. The 
core requirement for the general case would be, 


Z12,000 +4 216a + 80b + 1584c = bytes of core required. 


For a sample calculation of the core requirements 
consider the example of a spherical cap with 40 stations 
along the meridian, and an asymmetric analysis with two 
harmonics. Therefore, 


a= 4Q0(stations) x 2(harmonics) = 80 
b= 80 + 2 x 2(harmonics) = 84 


2 (harmonics) 


Q 
ul 


Thus, for the variables P, DEE, DST, 
3 x (4x4x80) = 3840 (words) x 4 = 15,360 bytes 


for the variable X, 
4 x (80) = 320 (words) x 4 = 1280 bytes 


for the varibles PHIXB, PHITB, 
2 x (80) = 160 (words) x 4 = 640 bytes 


for the varibles Z, ZO, 22, 23, ZDOT, 
5 x (4x84) = 1680 (words) x 4 = 6720 bytes 


mistiy, £0r the waribles ID, JD, IS, JS, 
G4 x (99x2) = 792 (words) x 4 = 3168 Lytes 


Therefore, the total size of the main dimension 
statement would be 27,168 bytes. This figure would be 
rounded up to the nearest even thousand bytes, ie. 28,000 
bytes. Finally, the core requirement for this example 
problem would be 


272,000 + 28,000 = 300,000 bytes. 
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III. IMPROVED POLE ROUTINE 


The SATANS code is based upon Sander's geometrically 
nonlinear equations under the conditions of small strains 
and moderately small rotations. The formulation is in four 
‘second order nonlinear partial differential equations in 


Teens Of UO, WV, W, and 4, where U, V, and W are the 
s 


meridional, circumferential and normal displacements 

respectively, and M is the meridional bending moment. The 
s 

nonlinear partial differential equations in the coordinates 


S, © , and t are reduced to uncoupled sets of linear 
differential equations in s and t by expanding the variables 
in trigonometric series in the circumferential coordinate 9Q, 
and treating the nonlinear terms as pseudo loads. The first 
and second derivatives in the meridional coordinate s are 
replaced by the conventional central finite difference 
approximations, ie. 


‘Zt = 1724s. { {2}. re a ey, (1) 
pi 1+1 sheng 
and 
x 2 
ato 17D (Cle). wwe iz}. + {Z} - |) (2) 
ae +1 7 oa ae 
th 
mere  f7) wS the vector of J, V,.W,.and Mat the.i 
a Ss 


station, & is the uniform dimension between stations, and 


primes denote partial derivatives with respect to “s. 


Applying these approximations to the governing set of domain 
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equations leads to 


Rea i2i jreoeizs TA zh: etsy: (3) 
1 1-1 1 z 5 1+1 a 


When the shell does not have a pole, fictitious stations 
one increment off of the shell are introduced at each end. 
Both the governing domain eguations and the boundary 
conditions are applied at the two boundary points. Thus, 
all finite difference approximations to the derivatives, 
including those of the boundary conditions, are of order 


Z 
we However, prior to the development of SATANS-ita, the 


treatment of the conditions to be applied at a poie at 
either end of a shell was handled by a simple Euier forward 
or backward difference approximation to the EaSse 
derivative, with truncation error of order 4. For example, 
for a pole at s= 0, where i= 1, the first derivative at the 


pole was approximated with 


= 1 - 3 4 
(z¥ Te AV = ty) (4) 


At the time this procedure for handling the pole 
conditions was developed (1967) it was thought that this 
would not significantly alter the solution. However, it has 


since been discovered that such is not the case. 


For the new pole routine, an expanded forward difference 


2 
approximation cf order A is used at s= 0 which takes into 
account the two stations after the pole, instead of just one 
station after the pole as in the Euler scheme. This 


approximation is 


ee= 1720) (302) 00+ 4 zy. = {Zz} _)- (5) 
1 1 2 3) 
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The conditions to be imposed upon the dependent 
variables at a pole are derived in Reference 14. They are: 


For N= 0, u =v =w =m = 0. 


Applying equation (5), these conditions can be put into the 


Batrix form 


10 0 0 u 000 0 u 0000 u 
0 420.0 v + 0co 640 Vv + 0077070 Vv = 0. 
0 0-3 0 W 004 0 Ww 0 0-1 0 W 
oO 0-3 n 000 4 m 0*0 0-1 n 
Sion SS 2 sy 3 


where the above 3 matrices are DL, DG, and DF within the 


SATANS programs. 
For N= 1, GV. =f =v = 2 = 0, 
where the plus sign applies at an initial pole, and the 


minus sign at a final pole. The matrix ferm for these 


conditions is 


=—3 0 0 0 u 400 0 u -1 000 Dt 
f1 0-0 v + 0 0 0:40 Vv + 00 0 0 v = 0. 
o 90. 1 0 Ww 0.0.0.0 W Q 0-00 
2.0.0 1 n 000 0 m 0 0 0.0 
Se See sg 
For N=2, u= v= w= m = 0 
s 
the matrix form is 
100 0 u 000 0 u 0000 u 
© 41.00 v + 000 0 v + 0 9 0-6 Vv = 0. 
O10 1.0 w O40"0 0 W oO: 0 0 W 
00 0-3 D 000 4 m 0 0 0-1 m 
Ss 1 si 2 s-i13 
Bor N > 2, u= v= w= = 0 
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and DL= identity matrix, DG= DF= null matrices. 


The sclution procedure in SATANS is an eiimination 
scheme and starts with 


ae * CEA Aza : ee ’ (6) 


where the values in ced based upon the Euler approximation 
are defined in Reference 14. The higher order approximation 


defines a new ie This new ae is obtained by 


Simultaneously solving the pole conditions 
[DL] cal + [DG] cae PCD] Te = AO Tie (7) 


and the domain eguation at station 2 next to the pole 


Cc B A = 7 8 
I a at Ces al ee {9}, (8) 
to eliminate let Thus, 
=f Aq = [fc - [B : 9 
{Zz} 3 { ys ( {93 3 [ He {Z} ; [ 1, {Zz} 1 (9) 


Substituting equation (9) into equation (7) gives 
+4 
DL DG DF A = iG = 
[ I z3e° sal gztes tik Ipbaass (Co ishs CEals (25 e 


B =p) 10] 
[ Les ae (19) 


Combining like coefficients of the {z} vector leads to 
DL 1-=[, DE, LA CC + DG) - [DF] [a] 'L[B 
EDL Lesh DE, b It 1) {Zhe Gide) sel DF dal iit })) 


== FOr TT Aa : 11) 
coe 5 jg Li, {3} { 


hinalLly, solving for oe yields 
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Zz == €pL— DE x A’ x Cc DG\= DF x A‘! xB Zz 
{ 7 [ aS [ oe { es 


Di = pF x Ae x. Cc - DF x a’! i 12 
AE Sd eal al EEN (12) 


-1 
Thus P = = [DL = DF x A* x Cc DG — DF x Av xB and 
Pan IF [ > oo if 3 of 


-1 
x = PDL =-opRiz°at xe - DF x aA’! <* The new 
{ }, [ : ou [ 51 {9}, 


[ PI, matrix has been placed into the "PMATRX" subroutine of 
SATANS-IIA and the new i, vector has been placed in the 


"FORCE" subroutine. 


A listing of the pole routine may be found in Appendix 
D. To anecEeperate thas new routine into a SATANS-I or-II 
program, first proceed to the "PMATRX" subroutine and remove 
the fifteen cards that are between, but not including, 
SEE (NN.GT.2) GO TO 90" and "11 CONTINUE". These cards are 
located after statement number "14" and just before 
statement number "11". Replace the cards removed by the 
@nes listed in Appendix D which read from" C IN PMATRX ™ to 
" 90 M3=MN". Then proceed to the "FORCE" subroutine and 
remove statement number "10". Replace statement number "10" 
with the nine cards listed in Appendix D which read from "C 
Bye erOrceE"s to ™DO 11 T= 1,4" +. Also place " COMMON 
JIBLS/IBCINL, IBCFNL " into the common area of the " FORCE " 


subroutine. 


This completes the implementation of the new pole 
routine into either SATANS-I or II. 
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IV. PROBLEM DESCRIPTION 


The geometry of the shallow spherical shell used in this 
study is identical to that used in Reference 1. Briefly, 
the shallow shell can be specified by the non-dimensional 
parameter X , where 

2 «atz4 if2 

Op 263 (1.-V ) | ( H/ h) : (1) 
Heats “the “rise “of the shell;ch is the thickness, and Y is 
Poisson's ratio. The mass density of the shell is a. All 


shells analyzed had the foliowing dimensions; 


Radii of Curvature R = Ro = 250 inches 
s 

Thickness h = 0.25 inches 

Modulus of Elasticity E = 30,000,000 psi 


Poisson's Ratio 053 


All buckling pressures obtained will be listed as a 
percent of the classical buckling pressure of a complete 
sphere, q4’ where 


2 eels 2 
47 PaZeL? feat Heres (CP TE- VEE) | (2) 
Forty stations were used over the meridian. The 
nondimensional time increment § t, where 
A 1/2 
C= It / < a 7%) ’ (3) 
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was taken as 0.05 for 3000 time steps, which is a total 
nondimensional time of 150. In addition, the axisymmetric 
analysis was repeated with a larger time step of 5 t= 0.2 
for a total time of 600. In this study m was selected such 
that t is equal to IT. The necessity for the long response 


time is explained in Reference 6. 


In the axisymmetric analysis only the N= 0 harmonic is 
considered. However, in the asymmetric analysis a_ second 
harmonic is excited by applying an incremental load in that 
harmonic. In addition, analyses of the shells X = 6, 7.5, 
and 11 were made using five harmonics. The step pressure 


load for the axisymmetric harmonic is 


[aren nS 2 q) {1}, (4) 
and the step pressure load for the asymmetric second 
harmonic is 

(n) (n) 


CS ee qo Ee lls (5) 


where n> QO, and a is taken as 0.0001. The value taken 
for the second harmonic in the asymmetric analysis was the 
Same as the critical harmonic for the static buckling 
analysis presented by Stilwell and Ball [2]. When there was 
an uncertainty as to which was the critical static harmonic 
the two harmonics in question were both tested. Run times 
using SATANS-IIA with a two-harmonic analysis for 3000 time 
steps and 40 stations on the meridian took an average of 28 
Minutes on the IBM 360/67. 


The parameter used to determine the minimum ioad at 


which dynamic buckling occurs is the peak value of V, called 
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Vv , where ¥ is defined as 
MAX 


"5 (0) G 
v= [ra acy J rf ar (6) 
0 t) 


ris the normal distance from the axis to the shell, a is 
B(O)e 
the maximum value ofr, W is the normal displacement of 


the axisymmetric response and § is the vertical distance 


from the base plane to the undeformed shell. The Vis a 
measure or the volume ort the shell deformation. The Fortran 


statements computing V and oe are given in Appendix E. 


When working a problem that requires these calculations the 


Nineteen cards are inserted directly into the "DYNAMIC 
subroutine right after the "IF" statement that calls tke 
"OUTPUT" subroutine. 


For convenience, the response in each asymmetric 


(0) 


harmonic is also measured uSing equation (6), with W 
(n) é 
replaced with W - The parameter V for the asymmetric 


harmonics does not represent a volume of deformation as it 


does for the axisymmetric harmonic. It can, however, be 
used to indicate the relative excitation of the asymmetric 


harmonics. 


The buckling criterion for both the axisymmetric and the 
asymmetric dynamic buckling analysis defines the critical 


load as that load P where a very small increase in P causes 


‘a very large increase in V - This is the same criterion 
MAX 
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Vo RESULTS AND DESCUSSION 


A. STATIC AXISYMMETRIC BUCKLING ANALYSIS 


Table I presents the new results from the static 
axisymmetric buckling analyses for A = 4 through 13 using 
the new pole routine. The two upper curves in Figure 1 
present a comparison of the new results obtained by 
SATANS-IIA with those obtained by Stilwell and Ball [2] 
using the SATANS-I program. AS can be seen in this figure, 
fairly significant changes in the buckling load occurred in 
the neighborhood of » = 4,5, and 9; and somewhat smaller 
differences cccurred in the region A = 10 through 13. The 
upper data points in Figure 2 present the comparison of the 
hew results from SATANS-IIA with those obtained by Huang 
[4]. This comparison shows a very good agreement between 
the two sets of results, except for the largest values of X. 
The new results have eliminated the differences that existed 
between the SATANS-{ results and Huang's results. 


B. DYNAMIC AXISYMMETRIC BUCKLING ANALYSIS 


Figure 3 presents the new results for the peak valve of 


v versus P for the various values of » tested. Table 
a (pt presents all of the new results for the dynamic 


axisymmetric buckling load. These loads are selected from 
figures constructed just like Figure 3. In every case, 
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except for *’°= W;°a' value of P slightly above the Bee 


value caused a V indicative of buckling, as well as a 
M 


nonconvergence of the iterative solution procedure. 


The lower two curves of Figure 1 present a comparison 
versus of the new axisymmetric dynamic buckling results 
with the previous buckling results obtained by Ball and Burt 
{1}. In every case the new critical pressure is lower than 
the critical pressure obtained using the Euler approximation 
at the pole. 


The lower data points of Figure 2 present a comparison 
of the new results with those obtained by Huang [5], by 
Seephens and Fulton [6], and by Stricklin [8]. Just as in 
the case of the static axisymmetric buckling analysis, the 
hew cesults compare much more favorably with the other 
results than did the results of Reference Ue. Pes 
interesting to note that the new results now tend to be 
Slightly lower than the other results, whereas the resuits 


of Reference 1 were aigher for almost all values of x. 
C. DYNAMIC ASYMMETRIC BUCKLING ANALYSIS 


Table III presents the new results for the critical 
pressures obtained from the dynamic asymmetric analysis. 
The second harmonics, or critical static harmonics, used in 
the analyses are also presented in Table IIi. A comparison 
of the critical pressures from the asymmetric analyses, 
Table III, with the critical pressures from the axisymmetric 
analyses, Table II, reveals that only the shell A =6 
buckled at a load below the axisymmetric buckling load. For 
the shell A = 7 the critical buckling load was slightiy 
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larger when asymmetric motion was considered. In all other 
cases the buckling was not influenced by the presence of the 
second harmonic. These new buckling results and those by 
Ball and Burt [2] are plotted in Figure 4. The new results 
can be seen to be significantly different from the SATANS-I 
results, where the asymmetric buckling loads were lower than 
the axisymmetric loads for five out of the ten values of 
tested. 


Except for » = 6 and 7, the relationship between aaa 


and P for the N= 0 harmonic, in the two-harmonic analyses, 


was found to be essentially identical to the relationship 
found in the axisymmetric buckling analysis shown in Figure 


=e Table IV A presents the v versus P data for both the 
MAX 


N= 0 harmonic and the second harmonic, for all values of A 


tested, except for XX = 6. Note that, except for X=7 and 


at, Vv for the asymmetric harmonic is generally very 
MAX 


small, even when the ¥ for the N= 0 harmonic indicates 
MAX 
that the shell has buckled. Thus, except for the shells A= 


6 and 7, the presence of the asymmetric motion does not 
influence the axisymmetric motion, and except for the shells 
» = 6, 7 and 11 the asymmetric motion is very small prior 


to buckling in the axisymmetric harmonic. 


A more detailed analysis of the shell X = 6 has been 
conducted since it was the only shell that revealed any 
Significant axisymmetric sensitivity to asymmetric motion. 
This shell was studied using two two-harmonic anaiyses (N= 
0, 1and N= 0, 2) and a five-harmonic analysis (N= 0, 1, 2, 
3, and 4). Figure 5 and Tables IV B and IV C contain values 


er Vv versus P for both of the asymmetric harmonics, N= 1 
MAX 
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and N= 2, in the two two-harmonic analyses, as well as the 
values of Re for the axisymmetric harmonic, N= 0. Figure 
6 and Table IV D present the values of y ae versus P for the 


N= 0,1,2,3, and 4 harmonics from the five-harmonic study. A 
comparison of the critical buckling load predicted from the 
results of the two two-harmonic analyses in Figure 5 with 
the critical load from the five-harmonic analysis obtained 
from Figure 6 shows that the presence of the additional 
harmonics results in the shell buckling at a slightly lower 
load (0.50), with significant motion in the N= 1 harmonic 
instead of the N= 2 harmonic (see the norconverged solution 
aor =  O.51), Which is the critical harmonic for static 
asymmetric buckling. Studies using five harmonics have also 
been conducted for A = 7.5 and A= 11. As can be seen in 
faple £VY D the critical harmonic for A= 7.5 remained N= 3; 
however, significant motion occurred in that harmonic at P= 
wigeand -44. In the case -of A = 11, relatively large 
asymmetric motion occurred in the asymmetric mode of N= 5 


Vice 6 at a value of P= .46. 


The comparison of the new results for the critical 
pressure for dynamic asymmetric buckling with those obtained 
analytically by Strick Lin Ee, by Akkas [9], and 
experimentally by Lock et al [7] is iliustrated in Figure 7. 
The comparison reveals an agreement with Stricklin in every 


case, in general a higher value of ace than those obtained 


by Akkas, and most importantly a very good agreement with 


Lock's experimental results. 


When making the comparison between the new results and 
those obtained by Akkas, it is necessary to look at the 
differences in the problem solution parameters used in the 
two studies. For example, buckling results obtained from 


SATANS-IIA using the same time increment as used by Akkas, 
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§ t= .2 for 3000 time steps, were significantly higher 
than those using the time step of &t= .05 for many values 
of A. Furthermore, the new resuits had , in some cases, 
instances of buckling occurring as far out in time as 130. 
Akkas, to shorten computer run times, observed the cap only 
for a time of less than 5. Furthermore, only the harmonics 
B="? or 2 or 3 were studied by Akkas for shells A= 5 
through 12. KE the critical harmonic is not studied, the 
predicted load will be too high. Thus, it appears that 


Akkas' lower bound loads may not be true lower bounds. 


Two additional features of the shell response should be 
noted. heest sheltis Al = 6G, 725, and 11 exhibited a 
non-buckled response in the axisymmetric harmonic to a load 
larger than the defined critical buckling load. This can be 
seen in Tables IV A and IV C. Second, and most importantly, 
the buckling load proposed by Ball and Burt [1], and used 
here, defines buckling to occur when the ae in the 


axisymmetric harmonic undergoes a large change due toa 


small change in P. Another criterion for dynamic buckling 
in the asymmetric analysis discussed in Reference 1 is to 
define the buckling load as that threshold load that 
initiates significant growth in the asymmetric harmonic. 


Re-examination of the V versus: P data an Table IV A 
MAX 
through D reveals that shells A= 6, 7, and 11 exhibited 


relatively large asymmetric motion at loads smaller that the 


defined buckling load when compared with other V values 
MAX 


for those shells, even though the numbers themselves were 


Small when compared with the axisymmetric harmonic. Shellsa 
<7.5 and 12 appear to be borderline cases. If the alternate 
criterion for buckling is used, the critical buckling loads 
Boeeshelis A= 6, 7, and 11 become 0.47, 0.45, and 0.45, 
respectively. The shells \ = 7.5 and 12 could have buckling 
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loads as low as 0.40 and 0.44, respectively. These values 
are more conservative than the definition based upon 
axisymmetric response. These five shells are the same five 
shells that exhibited an asymmetric buckling load lower than 
the axisymmetric buckling load in Reference 1. 
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VI. SUMMARY AND CONCLUSIONS 


A digital computer program for the geometrically 
nonlinear. analysis of totally arbitrarily loaded shells of 
revolution (SATANS-II) was modified to more accurately 
account for the conditions at the pole of the shell. This 
program, called SATANS-IIA, was used to determine the 
buckling load of shallow spherical shells of various sizes 
when subjected to static axisymmetric, dynamic axisymmetric, 
and dynamic nearly axisymmetric step-pressure loads of 
infinite duration. The cap sizes ranged from A= 4 to 13 
including \ = 7.5. A comparison was made between the new 
buckling results with the improved pole handling routine and 
the results that did not have the new pole routine. The 
comparison revealed a Significant change in buckling 
pressures, due solely to the change from an order A _ finite 
difference approximation of the first derivatives at tne 


2 
pole to an approximation of order A . These new critical 


pressures are in very good agreement with the results from 


other studies of the same spherical shells. This good 
agreement with other results, which came about as a result 
of the modification of the pole handling routine, is a 
strong indication that the manner in which the pole 
condition is handled is vital to the accuracy of the 


solutions obtained. 


In the asymmetric analysis, two harmonics were included 
for most of the shells; the axisymmetric harmonic and one 
asymmetric harmonic. Five-harmonic analyses were conducted 


for three of the _ shells. Two buckling criteria for the 
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asymmetric aralysis were considered. One defined buckling 
as that threshold load that caused a large increase ina 


deformation farameter, Me ae in the axisymmetric harmonic. 


The other, more conservative than the first, defined 
buckling as that threshold load that caused a large increase 


in the ae value for the asymmetric harmonic. Both values 


have been presented. 


The new static axisymmetric, dynamic axisymmetric, and 
even the dynamic asymmetric critical buckling pressure loads 
appear to te fairly reliable results for perfect, shallow 
shells. The effect of realistic imperfections remains to be 


determined. 
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1. TABLE IV A. Two-harmonic analyses for all values of 
A exept A= 
mm = 5 N= 0 and 2 N= 0 and 1 


ES] eae | ee aa | ae a 
i= 0 |. 
ju= 2 | 


Ace aaa 

[1659 | «1676 | .6606 [| w= 0 | . 1675 [6606 | 

0004787 | .o000566 | .0687 || x= 1 |.0003145]. 0687 | 
[aie are ae aa 


46 
ee ee ee ee 
= 1.001092 
=a, N=-0 .and.3 


a ec a 
of sess [poner [over [soos [sone 10s [1055 
= 3 

[ass fass [oss 
peo ne | nae [ ar 


= 5, N= 0 and 3 


as es 
IN= 0}.0703 .07228| .07429 2636] 07837 . 2076 
= 3].0001094|.001296|.0004304 | .002754 | .001188 | .000338 


ae 
m=O] .200 | 


40 


a 


lev Ife 2% “ern paaonaad~ovt 4 ; 
T Saa- 0 = : 

- ea). Bike T | 4, | - £h 
eee ee ee eee 
bo.) 2a7 ee Aad | 2020000, [NSEegy 

| u. 
Seorod. i! 
— é fae 0 me 
eee eee es ee 
Soe] ioe. ieee oa 
ee as 4b 108806 
ee BT Ps | 
3 MAS «: | eSar. BOSTG, ais 
| ae 


-05893 | .0607 .1713 |.1957 
. 0000566 | .0000703 | .0000364 |. 0000333 |.0000274 |. 0000299 


N= 0] .04738 |.04875 . 1576 i neon 
w= 4|.00003597|.00004635|.00004497 || N= 5 |.00008385 
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2. TABLE IV B. Two-harmonic analyses with N= 0 and 1, 
A-="6~onLys 


peumoiIseeni|ia ens). spuocsous | posva7 [sz 7) 53 
fu= of .4350 | .1385 | .1421_|-vweo | 1009 | 0053 


3. TABLE IV C. Two-harmonic analyses with N= 0 and 2, 


[n= 2.000234 0107 |.0293 | .01976 |.0a760|. 1223 | . 1419 


4. TABLE IV D fFive-harmonic analyses for selected 


shells. 


m= 6 N= 0,1,2,3, and 4 


|e ECE UE TS 
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= 2[.0187 |.0953 [.1059 [.1507 |.1502_|.1279 [2305 
54 

n= o|-3050 

= 1].3953_ 
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